function xdeltastCOMFn = deltastCOMFn(x)
xdeltastCOMFn = (-24730236602*x(1) - 11971114479*x(2) + 2094788099*x(3) - 457593477*x(4))/ ...
 24730236602 ...
;